Van der Waals interaction between two crossed carbon nanotubes 
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The analytical expressions for the van der Waals potential energy and force between two crossed 
carbon nanotubes are presented. The Lennard- Jones potential for two carbon atoms and the method 
of the smeared out approximation suggested by L.A. Girifalco were used. The exact formula is 
expressed in terms of rational and elliptical functions. The potential and force for carbon nanotubes 
were calculated. The uniform potential curves for single- and multi- wall nanotubes were plotted. 
The equilibrium distance, maximal attractive force, and potential energy have been estimated. 
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II. ANALYTICAL APPROACH 



A. Model 
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The van der Waals (VDW) interaction plays very im- 
portant role in Nano Electro Mechanical systems and 
Nano Electronic Devices^. Carbon nanotubes (CNTs) 
are promising material for creation a nanotweezers^ 4 -, 
nanoswitches^&i, bearings^^, nanotube random ac- 
cess memor y 11 : 12 , etc. The VDW forces are very criti- 
cal for understanding the growth mechanism of fullcrcncs 
and nanotube s 13 ' 14 : 15 and formation process of ropes 
and bundle o 16 ' 17 : 18 . Potentials for graphite layer s 19 ' 20 , 
two fullerene o 21 : 22 ' 23 : 24 , fullerene and surfac o 25 : 26 , nan- 
otube and surfac o 9 ' 27 : 28 , fullerenes inside and outside 
of nanotube o 29 ' 30 : 31 ' 32 : 33 are well studied. There are 
a number of publications devoted to interaction be- 
tween the inner and the outer parallel tubes such as 
single- (SWNTs)^^^^, double- 3 !^, and multi-wall 
nanotubes (MWNTs^^il. 



The present work is dedicated to interplay between two 
CNTs crossed at arbitrary angle. Problems related with 
application of molecular dynamics and density functional 
theory for such kind of calculations are discussed for ex- 
ample in Refs. |H and HE 



To evaluate potential between two crossed SWNTs 
or MWNTs we apply the continuum Lenard- Jones (LJ) 
model suggested by L.A. Girifalco^ 1 -. 



The model potentials for the VDW interaction are 
based on empirical functions whose parameters are ob- 
tained from experiment. It is remarkable that they have 
been so successful in providing a unified, consistent de- 
scription of the properties that depend on the weak in- 
teractions between and among graphene sheets, fullerene 
molecules, and nanotubes. 



The LJ potential for two carbon atoms in graphene- 
graphene structure is 
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where r is a distance, A = 15.2 eV • A 6 and B = 
24100 eV- A 12 are the attractive and repulsive constants 
respectively^ 9 -. We approximate the potential between 
two crossed SWNTs by integration of LJ potential 



v 2 J tp(r) dSi (f£2- 
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Mean surface density of carbon atoms for hexagonal 
structure is 
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0.393 atom per A , 



(3) 



where ai = 1.42 A is observed value of the C-C bond 
lengths for periodic graphite^ 7 -. If we know van der Waals 
interaction between two SWNTs then we may obtain in- 
teraction between MWNTs by summation over all pairs 
of layers 4 ^. 



B. Ancillary integrals 

To calculate the VDW interaction between two nan- 
otubes we have to take few useful integrals. For integral 
of LJ potential between two straight lines 



Iu{r) = / ip(r) dh dl 2 



we obtain 
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FIG. 1: Schematic drawing of interaction between: a) SWNT 
and line; b) two MWNTs 



where 7 is angle, r is distance between lines, and index 
II means "line-line" . For next integral between line and 
tube we have 

Iu(r,t)=t [ I u (r -t smfldp, r > t, (6) 

J — 7T 

where t is tube radius, r is distance between line and 
axis of tube, and index It means "line-tube". Figure 
la illustrates the schematic image of SWNT and line. 
Introducing new variable u — tan(/3/2) and using method 
of partial fractions we get 



Ilt(r,t) 
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where 



_ r 7^(2^+3) 

k — - , G A (x) = 7 — , 

t (x 2 -l)i 
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and using (7) we have 

tptt{r,ti,t 2 ) = v h / ===== rix. (10) 



(11) 



After some transformations we write 



I A - g A B ■ g B 



+ 



sin 7 \ r 
where multipliers for attractive and repulsive terms are 
n , s h l f bl+k 2kG A (>c) 

, Ay jrjw dx . (13) 



5 J bl -k y/k 2 - {x-bif 

The obtained Eqs. (12, 13) represents the elliptic inte- 
grals. In modern mathematics elliptic integral is defined 
as integral J R{x, y)dx, where R(x, y) is rational function 
of x and y, and y 2 is a qubic or quartic polynomial in x. 
With the apropriate reduction formula every elliptical in- 
tegral can be expressed in terms of elementary functions 
and canonical elliptic integrals of first, second and third 
kind. The method of integration is quite complicated but 
well know n 42 ' 43 ! 44 ' 45 . We would like to present only final 
answer. For attractive part we have the dimensionless 
parameter 



g A = g A KK(h) + g AE E(h), 



where 
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C. Tube-tube interaction potential 

Figure lb schematically shows the interaction between 
two MWNTs crossed at right angle. Parameter a 2 = 
3.44A is the average distance between two layers in 
MWNTs^l, d is the gap between tubes. 

In particular case these tubes may consist only of one 
layer. We notice t\ as a radius of first SWNT, and t 2 as 
a radius of second one (r = d + t\ + t 2 ). The interaction 
potential between two SWNTs is 



<ptt(r,ti,t 2 ) = 2v 2 



r+t 2 



Ilt(x,t 1 Wl + y' 2 (x)dx, (8) 



where y(x) = \Jt\ — {x — r) 2 , index tt means "tube- 
tube". 

Introducing dimensionless parameters 



x r r t 2 
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i,3=X..3 
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x (6162 + 62-61)*!, (16) 
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and matrixes of integer coefficients are 
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Analogically for repulsive part we write the dimensionless 
parameter 



9b = 9bkK(1i) + g BE E(h), 



where 
9b k = 



2 



(18) 
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The dimensionless coefficients for attractive and repulsive 
part of force are expressed as 



f2(i-l)j2(j-l) 
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Matrixes of integer coefficients {pbk} and {pbe} are 
placed in Table I and II respectively. 

As we see the final result (fTTj) is quite huge but it is 
working much better then usual numerical integration, 
because this analytical formula provides high accuracy 
and high speed of calculations. 



D. Tube-tube force 

The VDW resulting force caused by VDW interaction 
energy is 



F(r) 



d<ptt(r) 
dr 



Using expressions^ 3 - 

dK(x) _ E(x) K(x) 
dx (1 — x 2 )x x 



and 
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it is possible to obtain the analytical formula for VDW 
force. After usual differentiation over r we have 



F(r,h,t 2 ) 
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where 



!a = f AK K(h) + f AE E(h), 
Sb = f BK K(h) + f BE E{h). 



(24) 
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where 
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Matrixes {qbk} and {qbe} are given in Table III and IV 
respectively. 



E. Potential and force between equivalent tubes 

The potential and force can be expressed essentially 
simpler in the case when radii of interacted tubes are 
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TABLE I: Matrix {pbk} 
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7748055 


-2248848 


-1029636 

















2758140 


-938280 


-36368640 


27011520 


-1356976 













-1437660 


-31109610 


114941568 


-34271280 


-35205136 


5012854 












52840620 


-77089824 


-172119744 


193452624 


-11269956 


-5001276 














360818280 


-161478032 


-132247830 


47389620 


1327256 
















295052744 


-50123640 


-24636352 


1062586 



Symmetry 54162000 -2009850 -786849 

2858226 112076 
18436 



TABLE II: Matrix {pbe} 
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28878780 


-7265496 
























18117540 


-26421780 


-43621536 


16744920 




















4242735 


-322255500 


735299208 


-256914792 


-3437690 
















-66510360 


827431080 


-863979648 


-1027666080 


746908112 


-43635920 












104031900 


-508237800 


-2282192304 


5778243552 


-1895641688 


-527458544 


77409220 














4923519552 


-4510407600 


-4710377296 


4534092720 


-365372760 


-54614168 
















12018698404 


-4283796976 


-2131197060 


690518216 


8213905 


















5458820144 


-809953328 


-276097028 


11278540 



Symmetry 660421350 -23545932 -6367700 

24951224 777760 
114064 



TABLE III: Matrix {q BK } 

























-3356640 
492775164 









6971832 
-318885336 
-1357432104 
10893073968 








-5380200 
28477575 
2273850600 
-5600776440 
-7839981360 
20757703300 
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2109744 
79474500 
-991483500 
-1148679840 
10812253920 
-9567436112 
-6144532400 
8420997520 







-317520 
-40861422 
10951500 
1506696660 
-2776597920 
-1915609240 
7173492336 
-3892023180 
-1012896160 
973066230 





-42840 
4651668 
65362332 
-192184440 
-240160680 
1286509520 
-1300087984 
-138383496 
951956720 
-459154660 
-17001900 
37187456 



15624 
755433 
-1248828 
-11592924 
44486040 
-50195242 
-21785104 
107776844 
-103517008 
33353705 
9346988 
-8899384 
1348304 
155552 - 



TABLE IV: Matrix {<7sb} 


























68216904 










-48259071 
-1059904440 
21368220468 









11375658 
1921604958 
-11093841192 
-31611088152 
131005821528 
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8202411 
-1024742250 
-3163636875 
45974714520 
-55358109540 
-71790032408 
128236620500 







-7426692 
84564396 
4609741500 
-12201112980 
-25761560400 
87874942896 
-50711821008 
-32263986000 
33152425140 






2159073 
104520780 
-901760874 
-3041479980 
16321695750 
-14217370384 
-17650761720 
34135947792 
-13793389735 
-3602984700 
2709914346 





-124362 
-24489990 
-139540884 
861398916 
-594533340 

-3128825524 
6921793360 

-4488366288 

-1350545098 
3110622970 

-1197711564 

-49166500 
77206912 



-35469 " 
-1553454 
4928091 
17567676 
-110403234 
204320812 
-103004658 
-196048368 
380618007 
-270248958 
63736407 
24929412 

-17323848 
2281392 
236192 
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FIG. 2: Potential energies for interaction between pairs of 
identical SWNTs 
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FIG. 3: Potential energies for interaction between pairs of 
MWNTs of equivalent size. MWNTs contain exactly 10 layers 



equal t\ = t<i = t, therefore b\ = b 2 = b = r/t. Then we 
have for potential 

J ' A-g* A B-g 



sin 7 



+ 



where 



and 



g* A = g* AK K{2/b)+g* AE E{2/b), 
g* B =g* BK K(2/b) + g* BE E(2/b), 



9ak = - 



2ir 2 (5b 2 - 4) 



9ae — 



3(6 2 -4) 2 ' 

27r 2 (32 - 20b 2 + 116 4 ) 
3(6 2 - 4) 3 
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FIG. 4: Uniform potential for CNTs with arbitrary sizes. Dot- 
ted line is the universal curve suggested by L.A. Girifalco et 
al. 



2- 
0- 

z 

u: -2- 

-4- 
-6- 




Radius SWNTs 
t=6.79A (10,10) 
t= 1 . 1 8 A (15,15) 
t=13.57 A (20,20) 
t=16.96 A (25,25) 



I 

4 

o 

d,A 



FIG. 5: van der Waals forces between two identical SWNTs 
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FIG. 6: van der Waals force between pairs of MWNTs of 
equivalent size. MWNTs contain 10 layers 
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TABLE V: Dependence of potential well \iptt \ ( e V) from num- 
ber of layers for MWNTs of equivalent radii 



Radius (A) 


5 10 15 


20 


25 


50 


18.57 19.05 N/A 


N/A 


N/A 


100 


37.31 38.46 38.84 


39.00 


39.08 


150 


56.06 57.85 58.48 


58.78 


58.95 


200 


74.81 77.25 78.11 


78.54 


78.79 



9*bk = -7T 2 (460% 14 + 560386 12 + 321132b 10 
- 4736326 s + 18859526 6 - 3867648k 4 + 45107206 2 

- 2293760)/(1575(6 2 - 4) 8 ), (36) 

9 be = tt 2 (71296 16 + 972206 14 + 7634896 12 
- 15334246 10 + 77909446 s - 217561606 6 + 38781 1846 4 
- 400998406 2 -I- 18350080)/(1575(6 2 - 4) 9 ). (37) 

By totally the same way in the case of t\ = ti for force 
we have 



F*(t,r) 



4tt 2 



Af* A , Bf% 



sin 7 \ t 2 i 8 J ' 



(38) 



/ -=3^4^l / ^UJ +/ ^U )) ' (39) 



f% K = 13b 6 - 626 4 + 646 2 - 96, 



(41) 



f* AE = _256 6 + 366 4 - 1766 2 + 192, (42) 



f* BK = 97226 18 + 1296506 16 + 5376416 14 

- 58040366 12 + 114189766 10 - 509231366 s 
+ 1182118406 6 - 1805486086 4 

+ 1629716486 2 - 66060288, (43) 



f* BE = -147626 18 - 2851746 16 - 26599516 14 
+ 30296366 12 - 276227526 10 + 953263366 s 
- 2266119686 6 + 3515566086 4 

- 3218145286 2 + 132120576. (44) 



III. RESULTS AND DISCUSSION 

We have studied the WDV interaction between two 
crossed CNTs by using the continuum LJ approxima- 
tion. The analytical integrations for the potential energy 
of interaction between two identical SWNTs are plotted 
in Fig. 2. We assume nanotubes are crossed at a right an- 
gle in all our following illustrations and calculations both 
for SWNTs and MWNTs. We use parameter d = r - 2t 
to characterize the distance between tubes. Based on the 
results illustrated in Fig. 2, it can be concluded that real 
gap between surfaces of interacting SWNTs in equilib- 
rium state is changed slightly in range d = 2. 92. .2. 93 A. 
This distance is practically independent from the angle of 
nanotube intersection and the diameter proportion. For 
comparison in Ref. [2l] the equilibrium gap between two 
fullerenes Cgo is given as 2.95 A. 

In the case of MWNTs interaction we apply another 
attractive (A = 18.6eV-A 6 ) and repulsive (B = 29040eV- 

A ) constants which reproduce the layer distance of 3.35 
A and the elastic constant C33 = 4.08 GPa of graphite 37 . 
We assume that each pair of layers interacts as SWNTs 
and use summation over all pairs. The potential energy 
for two MWNTs of equivalent radii is plotted in Fig. 3. 
In these calculations we assume that each MWNT con- 
sists exactly of 10 walls. The equilibrium distance be- 
tween their surfaces is found to be do — 2.87 A, which is 
smaller then the equilibrium SWNT-SWNT gap. From 
our calculations follows that only several outer shells of 
MWNTs play essential role in the VDW interaction. For 
example, if two equal MWNTs with d = 200 A contain 
5, 10, 15, 20 or 25 layers then the minimum energy is 
-74.8, -77.2, -78.1, -78.5 or -78.8 eV respectively. Depen- 
dence of the minimum potential energy from number of 
inner layers for pairs of equivalent MWNTs is presented 
in Table V. 

It was found that the VDW interaction between Cqq — 
C60; Ceo-SWNT, Cgo-graphene, graphene-graphene, par- 
allel SWNT-SWNT, and parallel MWNT- MWNT can be 
described by a universal curv o 29 ^ 34 ' 40 . In our case uni- 
versal curve means that a plot of7p tt = tpu/\(Ptt \ against 
d = d/da gives the same curve for all tube-tube inter- 
actions, where ip^ t is the minimum energy and do is the 
equilibrium spacing for the two crossed tubes. As pointed 
above, the equilibrium distances are approximately con- 
stants both for SWNTs and MWNTs. 

We have calculated the minimum energy ip^ t for 
SWNTs of different radii (Table VI) as well as for 
MWNTs (Table VII). These results can be described by 
approximating formula 



CNT 



rtCNT 
°2 



t 2 



where C? WNT = -0.19285eV • A~\ Cl WNT 

-0.05847eV, &udCf IWNT = -0.388eV-A _1 
0.186eV are parameters for SWNTs and MWNTs respec- 
tively. It can be figured out from Tables VI and VII that 



(jMWNT _ 
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TABLE VI: Calculated depth |<p? t | (eV) for SWNTs (upper-right side) and approximation (lower- left side) 



Tube type 


(10,10) 


(15,15) 


(20,20) 


(25,25) 


Radius (A) 


6.79 


10.18 


13.57 


16.96 


(10,10) 6.79 


1.426 \ 1.434 


1.728 


1.979 


2.202 


(15,15) 10.18 


1.723 


2.080 \ 2.081 


2.384 


2.653 


(20,20) 13.57 


1.975 


2.385 


2.734 \ 2.731 


3.039 


(25,25) 16.96 


2.199 


2.655 


3.043 


3.388 \ 3.382 



TABLE VII: Calculated depth \ip%\ (eV) for MWNTs consisting of 10 shells (upper-right side) and approximation (lower-left 
side) 



Radius (A) 50 100 150 200 

50 19.03 \ 19.05 27.07 33.20 38.36 

100 27.04 38.43 \ 38.46 47.17 54.51 

150 33.17 47.14 57.83 \ 57.85 66.85 

200 38.33 54.48 66.83 77.23 \ 77.25 



approximating formula gives very good accuracy. Us- 
ing dimensionless potential Tp tt we can fit the potential 
of interaction between pairs of different SWNTs to one 
uniform curve and between pairs of different MWNTs to 
another one (Fig. 4). 

It is remarkable that plots for CNTs of the different 
radii fall down in the corresponding curves with accuracy 
of line thickness. For comparison Fig. 4 also shows a 
universal potential suggested by Girifalco L.A. et al.— . 

Figs. 5 and 6 show forces for two CNTs of equivalent 
radii. As we see in figures the behavior both for SWNTs 
and for MWNTs is qualitatively similar. The distance 
where attractive force reaches its maximum is in range 
3.40-3.41 for SWNTs and it is practically constant, 3.36 
A for MWNTs. 



proximation suggested by L.A. Girifalco to calculate in- 
teraction between two crossed CNTs of uniform and dif- 
ferent diameters. The exact formulas for potential en- 
ergy and van der Waals forces are expressed in terms 
of rational and elliptical functions. These formulas be- 
come essentially simpler in the case of interaction be- 
tween equivalent tubes. We estimated the equilibrium 
distance, maximal attractive force and potential energy 
for SWNTs and MWNTs. We plotted uniform potential 
for SWNTs and MWNTs. 
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